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f^ ' Abstract. Due to their small radii, M-dwarfs are very promising targets to search for transiting 

f^ , super-Earths, with a planet of 2 Earth radii orbiting an M5 dwarf in the habitable zone giving 

OQ ■ rise to a 0.5% photometric signal, with a period of two weeks. This can be detected from the 

1— I ' ground using modest-aperture telescopes by targeting samples of nearby M-dwarfs. Such planets 

3 would be very amenable to follow-up studies due to the brightness of the parent stars, and the 

' ' favourable planet-star flux ratio. MEarth is such a transit survey of ~ 2000 nearby M-dwarfs. 

OO ' Since the targets are distributed over the entire (Northern) sky, it is necessary to observe them 

individually, which will be done by using 8 independent 0.4m robotic telescopes, two of which 

(-H ' have been in operation since December 2007 at the Fred Lawrence Whipple Observatory (FLWO) 

C*) , , located on Mount Hopkins, Arizona. We discuss the survey design and hardware, and report on 

I I the current status of the survey, and preliminary results obtained from the commissioning data. 

O ■ 
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^—1 ■ 1. Introduction 

> 

vQ ' M-dwarfs are very favourable targets to search for transiting super-Earth exoplanets 

^— I ' in the habitable zones of their parent stars (e.g. Charbonneau & Deming 2007). This is 
CO ' predominantly a result of the very low luminosities of M-dwarfs relative to solar type 
^~~^ ' stars, which means the habitable zones reside at much smaller orbital distances. For 
t — the present discussion, we assume an M5-dwarf star, with mass and radius ^ 0.25 Mq 

^^ and 0.25 R0 respectively, orbited by a planet of 7 M^ and 2 R^, representing a typical 

target of the MEarth survey. Given these parameters, the M-dwarf has a luminosity 
~ 1/200 that of the Sun, and consequently, a planet receiving the same stellar insolation 
K>" ■ as the Earth would lie at 0.074 AU, corresponding to an orbital period of 14.8 days. The 
k>( I geometric transit probability for such a planet is also raised by a factor of approximately 
5_j ■ three, from 0.5% (for the Earth-Sun system) to 1.6 %. 

Cd ' The small radii of M-dwarfs also give rise to substantially deeper transits (0.5 % for the 

M5-dwarf, compared to 0.03 % for a solar- type host), and the combination of a smaller 
stellar mass and shorter orbital period increases the radial velocity semiamplitude, to 
~ 10 m s~^ compared to 1.3 m s~^ for the Earth-Sun system. 

The combination of these factors means that the detection of such planets around M- 
dwarfs by the transit method is feasible from the ground, using present-day observational 
techniques and detectors. MEarth is such a survey using modest commercially-available 
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equipment, based on a design study published in Nutzman & Charbonneau (2008) 



In this contribution, we summarise the results of the design study, present the novel 
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mode of operation of the MEarth survey, and describe the implementation and current 
status of the hardware. 



2. Survey design 

In order to detect a few super-Earth planets, or to place meaningful constraints on their 
incidence, for a geometric transit probability of ~ 1% we must survey a few thousand 
M-dwarfs. These should be bright in order to allow follow-up studies to be performed, 
and given the intrinsic faintness of M-dwarfs this translates to being very nearby: for 
example, an M5V star has an absolute J-band magnitude of Mj ~ 9. A reasonable limit 
for follow-up studies is J ~ 12, implying distances < 40 pc. 

Such nearby stars are expected to have high proper motions. We therefore appeal to the 
recently-completed LSPM-North catalogue (Lepine & Shara 2005), a survey of the entire 
Northern hemisphere based on photographic plates, which should be nearly complete for 
stars with proper motion > 0.15"/yr, containing a total of ~ 62 000 sources. We further 
restrict this to a sub-sample of 4131 stars within 33 pc (Lepine 2005) to remove high 
proper motion sub-dwarf contaminants from the sample. 

In order to select the M-dwarfs from the remaining population of (predominantly) main 
sequence dwarfs, we apply colour cuts using combinations of the T^-band magnitudes from 



Lepine & Shara (2005) and 2MASS J, H and Ks magnitudes, which are available for 
every star in our sample from the 2MASS all-sky data release (Skrutskie et al. 2006). 
The resulting "culled" sample consists of ^ 3300 probable M-dwarfs. 

These targets are spread in a relatively uniform fashion over the entire celestial North- 
ern hemisphere. We have therefore opted for a strategy of observing them individually. 
This mode of operation brings a number of benefits. In particular, the requirement on the 
field-of-view of the detector is substantially relaxed, since it need only be large enough to 
obtain sufficient comparison stars to perform differential ensemble photometry. It is also 
possible to tailor the parameters of each observation to the individual target M-dwarf, 
for example varying the exposure time to achieve the required signal to noise to detect 
planets of a given size, thereby saving observing time. 



Nutzman fc Charbonneau (2008), found that the most favourable targets for such a 



transit survey are, in fact, the smallest stars: although these are intrinsically fainter, 
the reduced count rates are compensated by having deeper transits, and their faintness 
increases the number of suitable comparison stars available for a given field-of-view. It is 
important to recall that for small field-of-view observations of single targets, the noise in 
the comparison light curve can become an important, or even dominant, contributor to 
the total noise budget. We therefore further choose to concentrate on the smallest stars, 
choosing 1976 with estimated radii < 0.33 R© (corresponding to a spectral type of ~ 
M3). 

Since the emission from M-dwarfs peaks in the near-IR, the number of detected stel- 
lar photons would be maximised by observing in this spectral region. Unfortunately, 
near-IR detectors remain prohibitively expensive, so MEarth uses conventional CCDs. 



Nutzman & Charbonneau (2008) found that the optimal passband for observing the (ex- 
tremely red) M-dwarf target stars was a filter cutting on longward of '^ 700 nm, and 
limited on the red end by the tail of the CCD quantum efficiency curve. The Ha line 
was deliberately omitted from the bandpass due to its potential variability in very active 
stars, which is obviously an undesirable feature for precision photometry. This passband 
approximates the sum of transmission of the Sloan i and z filters (Fukugita et al. 1996), 
and has been termed an "i -f z" filter. In reality, this is formed using a single piece of 
Schott RG715 glass, and is therefore economical and straightforward to manufacture, and 
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the final bandpass using an e2v CCD42-40 is very similar to the conventional Cousins / 
filter. 

Given these parameters, we can calculate the telescope aperture and field-of-view re- 
quired to survey our target stars. This was done assuming a standard noise model, includ- 
ing contributions from Poisson noise in the stellar counts and sky background, dark cur- 



rent, readout noise, and atmospheric scintillation using the formulation of Young (1967) 



Figure[T]indicates that a 40 cm aperture telescope with a 25' x 25' field-of-view is sufficient 
for the vast majority of the M-dwarfs with R < 0.33 Rq. 
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Figure 1. Necessary aperture vs. necessary field of view for LSPM M-dwarfs. Stars with radius 
> 0.33 R© are represented by open circles, while stars with radius < 0.33 R© are represented by 
filled circles. In order to determine the required telescope aperture, we demanded the photometric 
precision necessary to achieve a 3cr detection of a transiting 2 R® planet in a 150s integration 
through the i + z filter. The range of the x and y axes match that of the stars with radii 
< 0.33 Rq, while 30% of the M-dwarfs with radii > 0.33 R© require more than a 100 cm 
aperture, and thus fall above the plot limits. 

We have opted for eight such telescopes, which leads to an expected 2 yr duration to 
complete the survey of the habitable zones of these 1976 M-dwarfs. The design study 
indicates that a yield of 2.6 habitable zone super-Earths would be predicted if the true 
occurrence of these planets was 10% around our targets, with larger and closer-in planets 
being easier to detect. A null result would limit the occurrence of > 2 R^ super-Earth 
planets in the habitable zones of late-M dwarfs to be < 17% at the 99% confidence level, 
a result that again becomes a stronger limit for closer-in planets. 



3. Mode of operation 

In order to improve the survey efficiency, and in particular, to increase the number 
of M-dwarfs that can be monitored at once by each telescope, we have adopted a novel 
detection strategy (illustrated in Figure [2]). Routine observations will be carried out at 
extremely low cadence, with ^ 2 visits to each target per transit timescale (we assume a 
mid-latitude transit, which leads to a duration 0.866 times that of an equatorial transit). 
The cadence is further limited to not being less than once every half-hour to assure 
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a reasonable number of observations of each target per night, and to catch shorter- 
period transits due to planets interior to the habitable zone. We intend to detect transits 
while they are still in progress via real-time analysis of the images as they are taken. 
This information can then be used to immediately direct follow-up resources (e.g. other 
MEarth telescopes) to confirm or reject the hypothesis that there is a transit in progress. 
In practice, this will probably proceed by immediately obtaining several more data- 
points, to combat the effects of noise and light curve systematics, and then following 
the remainder of the event, if it turns out to be real, at high-cadence to obtain a well- 
sampled transit egress. By continuously re-evaluating the probability of there being a 
transit in-progress upon obtaining additional data-points, we can quickly reject false 
positives without any major effect on the remainder of the targets. We can therefore 
cope with a relatively high false alarm rate. 




Time (days) 

Figure 2. Illustration of a possible sampling strategy for transit detection with MEarth. The 
solid line shows a transit model. The dashed vertical lines indicate the sample times during the 
normal (low-cadence) operation up to the first point which deviates by > la in the downward 
direction. Immediate follow-up is triggered at this point, confirming the downward deviation. 
Sampling continues at high-cadence until after egress, at which point the sampling returns to 
normal until another significant downward light curve deviation is detected. 



4. FoUowup and false positives 

As with any transit survey, in order to confirm the planetary nature of any detections, 
radial velocity follow-up observations will be required. This should be possible with ex- 
isting instrumentation, with several groups demonstrating the required radial velocity 
precision on M-dwarfs using conventional instruments in relatively blue passbands. Ide- 
ally, in order to take advantage of the greater continuum flux of M-dwarfs in the far 
red and infrared, radial velocities would be derived from lines in these spectral regions. 
Unfortunately the conventional Iodine cell technique is not usable in the far-red, and 
the HARPS (and by extension, HARPS-NEF) spectrographs cannot be used beyond 
^ 690 nm. The new TRES spectrograph on the Tillinghast 1.5 m telescope at FLWO 
appears promising in this respect, since it can operate around the atmospheric window 
close to 900 nm, as do several IR radial velocity methods currently under development, 
for example Blake et al. (2007, 2008) and the T-EDI project (Lloyd et al., this volume). 

It should be noted that the target selection employed for MEarth largely eliminates the 
most common sources of astrophysical false positives that plague conventional wide-field 
transit surveys. In particular, the selection by proper motion eliminates giant host stars 
from the sample, and hierarchical triple systems are extremely unlikely due to the very 
late spectral types that they must have in order to be included in the MEarth sample 
(the brightest star in these would be a mid-M dwarf due to the colour selection we apply) . 
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Triples may also be resolvable with high-resolution imaging due to their close proximity 
to the Earth. Finally, blends with background binaries are substantially reduced, not 
least because our pixel scale of 0.75"/pix is ^ 1/20 that of a typical wide-field transit 
survey, and the high proper motions of our targets mean that previous epochs of imaging 
can be used to resolve many of these systems, by looking for the background source at 
the present position of the M-dwarf. 



5. Status 

At the time of writing, MEarth has 5 telescopes in operation, with two operating 
routinely since January 2008, and the final three since June 2008. We anticipate having 
the remaining three telescopes operational by September 2008. All of the telescopes are 
housed in the same building. We note that the use of multiple sites at different longitudes 
would be beneficial to improve phase coverage, but would be much more expensive, 
since we were able to re-use an existing building and much of the associated support 
infrastructure at FLWO, affording a very substantial cost saving. 

Data are presently reduced in real-time, although the real-time analysis software is not 
yet in-place, using a modified version of the pipeline from the Monitor project, described 
in Irwin et al. (2007) We have already achieved photometry at the ~ 0.6% level on 
M-dwarf targets over ~ 90 nights, and are currently investigating methods to further 
reduce the level of systematics present in the data. This has been sufficient to detect low- 
amplitude (few percent peak-to-peak) rotational modulations in several of our targets 
for which sufficient observations are available. 
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Figure 3. Example of an M-dwarf light curve from MEarth data taken during the first season 
of operations, plotted in / magnitude as a function of sample number (numbering from 1 for the 
first data point). Vertical solid lines denote boundaries between different nights of observations, 
which are numbered at the top of the diagram. Data were taken on 58 separate nights for this 
object. The estimated spectral type of the star is M4V, with an I-band magnitude of ^ 9.4. 
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Figure 4. Example of an M-dwarf light curve showing periodic photometric variations, presum- 
ably due to rotation. The upper panel shows the light curve plotted as a function of heliocentric 
Julian day number, and the lower panel shows the curve folded on the best-fitting rotation 
period of 11.661 days, plotted versus orbital phase, where 1 unit represents a single rotational 
period. The estimated spectral type of the star is M4.5V, with an I-band magnitude of ~ 10.5. 



A new M-dwarf eclipsing binary system was also discovered during the first few weeks 
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of routine MEarth observations. The system is a near equal mass pair, with estimated 
component masses ^ 0.3 Mq from photometry, and an extremely short orbital period of 
0.77 days. The light curve of this object obtained using one MEarth telescope is shown 
in Figure [51 




Phase 
Figure 5. Phase-folded light curve of the M-dwarf eclipsing binary discovered in MEarth data. 
The upper panel shows the /-band magnitude, over-plotted with an eclipsing binary model 
(solid line) generated using a modified version of jktebop (Southworth et al. 2004a,b; Popper 
& Etzel 1981; Etzel 1981) including out-of-eclipse modulations (presumably due to starspots) 
synchronised with the binary orbital period of 0.77 days. The lower panel shows the residuals 
(data — model) of the fit. Note the differing scales on the vertical axes in the two panels. 

Such eclipsing binaries are of vital importance for constraining the mass-radius relation 
on the main sequence, and hence theoretical models of low-mass stellar evolution, which 
are remarkably poorly-constrained even at old, main sequence ages for stars below 1 Mq. 
Indeed, many of the observed and well-characterised systems show substantial discrepan- 
cies with the predictions of the theoretical models (e.g. Ribas 2006). These uncertainties 
have a direct impact on the reliability of stellar parameter estimates, and hence on planet 
properties inferred from analysis of transiting systems. Indeed, as detailed elsewhere in 
this volume, many of the transiting planet observations have reached sufficient quality 
that this uncertainty in the parameters of the host star is the limiting factor in our knowl- 
edge of the properties of the planets. Such M-dwarf eclipsing binaries are straightforward 
to detect from MEarth data, and it should be possible to characterise them extremely 
well due to their brightness, allowing model-independent mass and radius estimates with 
uncertainties of only a few percent to be made (e.g. Andersen 1991). 
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